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^ , ABSTRACT 

! Morphologically it appears as if the Vela X PWN consists of two emission 

^ ! regions: whereas X-ray (~ 1 keV) and very high energy (VHE) H.E.S.S. 7-ray 

\C> . observations appear to define a cocoon type shape south of the pulsar, radio 

^ ! observations reveal an extended area of size 2° x 3° (including the cocoon area), 

^ ! also south of the Vela pulsar. Since no wide field of view (FoV) observations 

^ \ of the synchrotron emission between radio and X-rays are available, we do not 

Q \ know how the lepton (e^) spectra of these two components connect and how 

^ \ the morphology changes with energy. Currently we find that two distinct lepton 

^ \ spectra describe the respective radio and X-ray/VHE 7-ray spectra, with a field 

^ \ strength of 5/iG self-consistently describing a radiation spectral break (or energy 

maximum) in the multi-TeV domain as observed by H.E.S.S. (if interpreted as 
IC radiation), while predicting the total hard X-ray flux above 20 keV (measured 
by the wide FoV INTEGRAL instrument) within a factor of two. If this same 
field strength is also representative of the radio structure (including filaments), 
the implied IC component corresponding to the highest radio frequencies should 
reveal a relatively bright high energy 7-ray structure and Fermi LAT should be 
able to resolve it. A higher field strength in the filaments would however imply 
fewer leptons in Vela X and hence a fainter Fermi LA T signal. 
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Subject headings: pulsars: individual (Vela X) — Radiation mechanisms: non- 
thermal — ISM: jets and outflows — ISM: supernova remnants — gamma rays: 
theory 



Introduction 



Th e Vela-Puppis reg ion was first identified as a bright radio region by lMcGee. Slee &: Stanley 



(Il955l ). iRishbeth I (Il958l ) compared Vela X, Y and Z (the radio counterparts of the Vela SNR) 
and remarked that Vela X is the most prominent radio feature, with the specific distinction 
of having an almost fiat radio spe ctrum. A total Vela X radio energy spectrum of the form 

measured by lAlvarez et al.l (120011 ) , w hich was bas e d on fiux measure- 
ments up to 8.4 GHz. However, between 8.4 and 30 GHz, iHales et al.l (120041 ) resolved the 
filamentary and PWN structures and found that the spectra are overall harder: Whereas the 
PWN spectrum is F^, cx z/O i^O-OS^ ^^j^g filamentary spectrum is of the form F^, cx //-0-28±o.09^ 



Following EGRET observations of the Vela region, Ide Jager et al.l (Il996al ) showed that 
the observation of inverse Compton (IC) scattered HE 7-rays from the electron component 
producing this radio component, should allow us to constrain the maximum radio frequency 
and associated field strength. This electron component should IC scatter three target photon 
components: the Cosmic Microwave Background Radiation (CMBR) field, the Galactic far 
infrar ed radiation (FIR) field as a result of reradiation of dust grains, and the local starlight 
field JPorter et al.lboOTh . 



Further insight into the nature of Vela X was revealed by the detection of a cocoon of X- 
ray emission south of the Vela pulsar by ROSAT (IMarkwardt fc Ogelmanlll995l ). The X-ray 
data (with limited energy extent up to 2.4 keV) was sh own to be consistent with a thermal 
origin, but recently iLaMassa. Slane fc de Jager I (120081 ) used XMM-Newton observations of 
the southern head of this cocoon to show that the addition of a non-thermal synchrotron 
component with photon index as steep as ~ 2.3 is favoured by the data. By rescal ing the 



H.E.S .S. fiux relative to the XMM-Newton FoV, a field strength of 5/iG was found. iMilne 



( 119951 ) also identified a bright radio filament running along this cocoon and models for such 
filaments require a strong thermal component, pulsar- generated magnetic fiux and relativistic 
particles (IReynolds Ill988l ). 



Conclusive proof for the existence of relativi stic particles in the X-r ay cocoon was re- 
vealed by H.E.S.S. observations of the Vela region (jAharonian et al. II2006I ). In fact, H.E.S.S. 
detected a spectral maximum from Vela X, which allows us to probe the low energy spectral 
component, as well as the energy at which the spectral turnover occurs. An overlay of the 
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H. E. S. S. image is shown in Fig. 8 of lde Jager &: Djannati-Atail (120081). showing t he overlap of 
the VHE 7-ray cocoon and the bright radi o filament referred to by iMilne ( IQQsi) . Whereas a 



hadro nic in terpretation was suggested by iHorns. Aharonian &: Santangeld (120061 ) , Ide Jager . 
( 2007 ) and de Jager k. Diannati-Ataj (2008h r eviewe d evidence in favour of a leptonic IC 



origin. Recently iLaMassa. Slane fc de Jager I (120081 ) also found that the thermal particle 
density at the head of the cocoon where bright VHE 7-ra y emission was found, is only 0.1 



cm ^ - a factor 6 lower than the lower limit employed by iHorns. Aharonian fc Santangelo 



(120061 ). This places more stringent constraints on the brightness of a hadronic signal from 
Vela X. 



Mangano et al.l (120051 ) derived an X-ray photon index of 1.50 ± 0.02 within a distance 
of 0.5 arcmin from the pulsar. This index represents the uncooled (by synchrotron losses) 
electron spectral index injected into the PWN, since a steepening is observed for distances 
greater than 0.5 arcmin. The corresponding uncooled post shocked electron spectral index 
of 2.0 is also consistent with the pre-break spectral index detected by H.E.S.S., if the signal 
from the latter has an inverse Compton origin. 



de Jager I (120071 ) has shown t hat Vela X most likely consists of a 2-component non 



thermal spectrum: From Fig. 2 of Ide Jager I (120071 ) the lepton spectrum corresponding to 
the VHE 7-ray spectrum only connects to the radio lepton spectrum if the field strength of 
the radio structure is > 200/iG, if we extrapolate the uncooled lepton index of 2.0 towards 
the radio domain. How ever, the fiel d strength in the radio filaments is more likely to be in 
the range 10 to 50 /iG jMilne \\l99^ . 



The observation of an X-ray photon index > 2.0 in the cocoon by lLaMassa. Slane fc de Jager 



( I2OO8I ) may be due to synchrotron cooling inside the X-ray band (below 7 keV) and to test 
this we need hard X-ray observations with a FoV covering the full size of Vela X. Only 
INTEGRAL (20 - 300 keV) has a large enough FoV to cover the extended Vela X (and by 
default cocoon) region. However, only point-like unpulsed emission centered on the pulsar 
(consistent with the 12 arcmin PSF) with a photon index of —1.86 ± 0.03 was detected by 
INTEGRAL (L. Kuiper, personal communication, 2007). Thus, synchrotron burnoff must 
have terminated most emission above 20 keV at distances > 12 arcmin from the pulsar to 
levels consistent with the se nsitivity of INTEGRAL . This spectrum is consistent with the 
OSSE spectrum reported by Ide Jager et al.l (Il996bl ). covering at least the total cocoon re- 
gion. The BeppoSAX PDS, with a FoV of 1.3 ° (FWHM), also obser ved the Vela pulsar above 
20 keV, covering most of the cocoon region (IMangano et al.l 120051 ). Although the observed 
photon index of the nebula is slightly steeper (2.00 ± 0.05), the energy flux is comparable to 
that of INTEGRAL. 



In this paper we construct a time dependent post shocked injection spectrum Q{E,t) 
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(with E the lepton energy). This spectrum consists of two components (as motivated a 
above): Qu^E^t) as source for the extended radio emission and Qx{E,t) as source for the 
X-ray and VHE 7-ray emission, confined mostly to the cocoon region. In this paper Fermi 
Ly4r observations should be able to test the spectral connection between Qx and Qji. 

These two distinct source spectra are normahsed to the time averaged spindown power, 
giving conversion efficiencies r^/j and rix of spindown power to leptons. The time dependent 
transport equation with synchrotron losses is then solved to produce the present time electron 
spectra dN^/dE and dNx/dE. The field strength required for the synchrotron loss term is 
obtained from H.E.S.S. observations, if we assume that the VHE 7-ray signal is the result of 
IC scattering. Since H.E.S.S. measured the spectral energy maximum, we find the magnetic 
field strength which reproduces the observed multi-TeV spectral break, so that dNx/dE 
is uniquely determined in the energy range 5 to 100 TeV. We then show that this same 
dNx/dE (with B derived from the VHE spectral break) also predicts the total hard X-ray 
fiux from the PWN measured by INTEGRAL correctly within a factor two. For ?7/j < 1 we 
also show that dN^/dE (with the abovementioned B) reproduces the total integrated radio 
fiux, which in turn predicts a high energy 7-ray component, which can be tested with Fermi 
LAT. 



2. The double source spectrum of Vela X 

In the discus sion below we will negl ect (to a first order) the effect of adiabatic cooling 



on the spectrum (jPacini fc Salvatil Il973l ) since the reverse shock (applicable to Vela X) is 
expected to have adiabatically heated the particles during the PWN crushing phase. Energy 
losses due to adiabatic cooling during the pre-crushing phase is then partially cancelled 
during the crushing phase (Ferreira & de Jager 2008, in preparation). 

Let Qji{E,t) and Qx{E,t) be the respective source spectra corresponding to the radio 
and X-ray synchrotron components, and L[t) be the time dependent spindown power, of 
which fractions ?7r and rjx are converted to these two components. Based on the discussion 
in Section 1, we model these source spectra as 

QR{E,t) = Qo,i?,(t)(^/^max,K)-^ (1) 
QxiE,t) = go,x(t)(^/^max,x)-', (2) 

for -Emin,/? < E < E^s_y, ji and -Emin.x < E < -E'max.x in the two respective equations. These 
upper and lower limits are currently weakly constrained by observations. The radio spectral 
index of 0.39 requires a particle index of p = 0.39 x 2 + 1 = 1.78, whereas an X-ray (VHE 
7-ray) lepton spectral index of 2 would reproduce the uncooled photon spectral index of 1.5. 
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By requiring the total energy in each respective component to be a fraction rjji < 1 and 
rjx < ^ (with rjji + rjx < 1) of the spindown power, we can show that (with the minimum 
electron energy for the radio component much smaller than the corresponding maximum) 



QxiE,t) 



n (F (a^LMl 

QR[E,t) ~ I 

VxL{t) 




(3) 
(4) 



In (-Emax.x/ -Emin,Xj 

Neglecting adiabatic losses (as discussed above), but including radiation losses with timescale 
ri.ad(-E') = E/{dE/dt)j.i^^, where [dE / dt)^.^ is the sum of the synchrotron and IC (Klein- 
Nishina effects included) energy loss rates, the transport equation can be solved to give 
the present day (age T) tot al particle spectrum for each component i = i?, X as (see e.g. 
Zhang. Chen k Fang I J2008h ) 

T 



dN,{E,T) 
dE 



Qi{E,t)exp 



T-t 



dt. 



(5) 



10 \ ^rad(^). 

This expression effectively allows leptons to be accumulated over time T in the nebula. The 
resulting spectrum dN/dE is however modulated (steepened) at the high energy tail by 
radiation losses. 

Since radiation losses are negligible for the radio component, the integral over the source 
term Qpt yie lds the change iri rotational kinetic energy (AKEj-ot = "~^^)/2) since birth, 
as shown by Ide Jager I (120081 ) . so that (for p < 2) 

dNn{E,T) _rfE{2-p)AKE 



dE 



max,R 



-E- 



(6) 



The term r/^ is the time averaged conversion efficiency. Ivan der Swaluw &: Wul (1200 ll ) found 
that a pulsar birth period of Pq = 40 to 50 ms is required to reproduce the observed ratio 
(~ 0.25) betwe en the PWN and SNR radius. The recently measured pulsar braking index 
n = 1.6 ± 0.1 (jPodson. Lewis &: McCuUochl 120071 ) implies an age between T = 11 and 15 
kyr given the abovementioned range for Pq. This range of values (or exact choice of braking 
index) do not change our conclusions significantly. 



3. The particle, synchrotron and IC spectra of the cocoon 

The lepton spectrum for the X-ray/VHE 7-ray emission region (mostly the cocoon 
region) derived from Eqn.([5]) (parameters in Table 1) is constrained by the following infor- 
mation: an injection (uncooled) lepton spectral index of 2.0 as discussed in section [1], the 
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IC no rmalisation of the VHE cocoon spectrum as provided by H.E.S.S. (lAharonian et al 



20061 ). the multi TeV break energy observed by H.E.S.S., spindown from Pq = 40 ms to the 
current period of 89 ms within time T (~ 11 to 15 kyr), an imphed time averaged conversion 
efficiency of rjx = 3 x 10~^, and a field strength of 5/iG to reproduce the observed VHE 
break energy. 

It is clear that this model IC spectrum agrees relatively well with H.E.S.S. measurements 
(Fig. 1 - thick solid line at very high 7-ray energies). The CMBR mostly contributes to this 
spectrum, whereas higher energy contributions from the galactic FIR and starlight photon 
fields are suppressed due to Klein Nishina effects. However, in the next section we will 
see that the latter two radiation fields contribute significantly to the high energy IC 7-ray 
counterpart of the radio spectrum, where the IC scattering is in the Thomson limit. 

The lepton spectrum of the cocoon, if self-consistent, should be able to predict the total 
X-ray spectrum from Vela X. From Fig. 1 it is clear that the model synchrotron spectrum 
shows a radiation break around 1 keV (~ 2 to 3 x 10^^ Hz) and predicts the 20 - 300 keV 
INTEGRAL fiux within a factor 2, although with a steeper spectrum. The latter is expected 
to harden if we allow the field strength to increase towards the pulsar in a more sophisticated 
model incorporating a spatial dependence as well. 



4. The particle, synchrotron and IC spectra of the total Vela X PWN 

To first order we assumed th e same field strength (5 x 10~^G ) of the entire Vela X 



structure as found for the cocoon ( LaMassa. Slane fc de Jager 20081 ). This may not be the 



case, but by comparing Fermi LAT with high frequency radio observations we should be able 
to measure the field strength in this extended structure if the high frequency (radio) cutoff or 
turnover can be independently determined. In this case synchrotron cooling is negligible and 
the steady state spectrum is given by Eqn IQ. To reproduce the absolute radio flux (for a 
particle index of p = 1.78) we had to assume efficiencies off]R = 0.35 and 0.45 corresponding 
to the assumed spectral cutoffs 10^° Hz and 10^^ Hz respectively. The leptonic parameters 
corresponding to the latter cutoff is listed in Table 1 under "radio" . From Fig. 1 it is clear 
that the maximum frequency should be > 10^° Hz. 

While the lack of high-frequency radio measurements prohibits a clear constraint on 
the cutoff of this component, Fermi LAT observations can constrain the maximum of the 
spectral energy distribution of leptons in the radio nebula. As an example we employ cutoffs 
at the abovementioned two frequencies to calculate the corresponding expected high energy 
7-ray spectra. These predictions are shown as dashed lines in Fig. 1 (marked "HE 7- 
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ray IC"). The hard lepton spectrum followed by an assumed cutoff at an electron energy 
Ema.x,R = 7max,ij"^c^ should produce broad 7-ray features at three energies: 'i-^l^a.^ r^'^i 
corresponding to the temperatures = 2.7K (CMBR), 25K (FIR) and ~ 4000K to ~ 8000K 
(starlight). 



Although these energy densities are comparable in the local neighbourhood (jPorter et al. 



20071 ). we investigate deviations at the location of Vela X: The IRAS- Vela shell (IVS) of dark 
clouds and cometary globules centered on {i^^ = 259.85°, b^^ = —8.25°) wit h a radius of 
(5.1° ± 0.2°) and annular RMS width of ~ 2° is located at a distance ~ 400 pc (ITestori et al 
20061 ). With a lOOyum IRAS ring flux of ~ lOMJy/sr we obtain a luminosity of ~ 10^® erg/s, 
and assuming a distance of r = 100 pc between the IVS center and Vela X, we obtain an 
additional FIR energy density of 0.002(100 pc/r)^ eV/cm^ at Vela X, which is well below 
our assumed local averaged energy density. Other potential target photon fields fro m the 



Vela molecular ridge and Gum Nebula are located directly behind the Vela pulsar (ISahu 



I992I ) so that the contribution from these local clouds would be severely suppressed by the 
head-on-tail effect of the IC process, since the (1 — cos^)^ term (with 6 small) in the energy 
loss rate term is much less than unity. 

Furthermore, fluctuations in the local distribution of stars can also produce deviations in 
the starlight energy density at Vela X. To investigate this we extracted the Hipparcos catalog 
of ~ 118,000 stars, of which ~ 73,000 have parallax measurements at the > 3a detection 
level. This subset was used to calculate the starlight energy density at any 3D field point 
within ~ 300 pc from E arth. Within ~ 20 pc from Earth we obtain a value of ?7rad = 0.24 
eV/cm^, consistent with iPorter et al.l (120071 ). whereas a value of ?7rad = 0.44 ± 0.12 eV/cm^ 
is obtained for t he Vela PWN, with the error reflecting the uncertainty in the distance to 
the Vela pulsar (IDodson et al.ll2003l ). This mean value (0.44 eV/cm^) was used in Fig. 1, 
but assuming a mean temperature of 7500 K to account for the stellar types contributing to 
this excess. 

Whereas the Fermi L^TwebpageSonly lists sensitivities for point sources, we scale the 
actual EGRET upper limits for the Vela X e xtended source (extracted from the unpulsed 
skymaps) obtained by Ide Jager et al.l (Il996al ) by a conservative factor of 30 downward to 
account for the improved Fermi L^IT sensitivity in the same energy bands for an integration 
time of 1 year. The EGRET integral flux upper limits and GL AST /Fermi LAT sensitivities 
are shown in Fig. 1. 



^ http: / /www-glast. slac.stanford.edu/software/IS /glast Jat_performance.htm 
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5. Discussion and Conclusion 

In this paper we have shown that a pulsar birth period of ~ 40 ms (to reproduce the 
observed ratio of PWN to SNR radii), spinning down within 11 to 15 kyr to the present 
period of 89 ms while converting rjx ~ 0.3% of its spindown power to leptons (e"*^) in the 
PWN with a spectral index of 2, would reproduce the observed H.E.S.S. flux in the Vela X 
cocoon via the IC process. A nebular averaged held strength of 5//G would also reproduce 
the observed TeV spectral break, Furthermore, this same spectrum predicts the BeppoSAX. 
PDS and INTEGRAL total fluxes within a factor of two. 

Assuming the same 5//G held strength in the much larger 2° x 3° Vela X radio nebula, 
we need to invoke a second e^ component with electron spectral index 1.78 and a much 
larger conversion efficiency rjR ~ 0.4 to reproduce the observed total radio spectrum up to 
at least 8.9 MHz and possibly 30 MHz. This spectrum should cut off or steepen signiflcantly 
to meet the less energetic cocoon spectrum. 

The implication of this second (radio) component is the prediction of an IC component 
at HE 7-rays covering the same 2° x 3° area south of the pulsar. The visibility of this 
component in the HE 7-ray domain for Fermi LAT is further improved due to the existence 
of the local FIR and starlight target photon fields, scattered by the highest energy leptons 
in this second component up to ~ 10 GeV 7-rays (in the Thomson limit) where Fermi LAT 
has good sensitivity. 

If the electrons radiating the radio component are trapped in compressed regions (e.g. 
filaments) with field strength ^ 5fiG, we would require fewer leptons and hence lower ef- 
ficiency rjji to reproduce the radio emission. A HE 7-ray fiux measurement, or upper limit 
thereof, will provide constraints on this field and hence rjR. 

PO Slane acknowledges support from NASA contract NAS8-39073. The referee is 
thanked for useful comments. 
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Table 1: Parameters associated with the two lepton components of Vela X for a pulsar birth 
period of 40 ms, and an assumed braking index of 2.5 to give an age of 11 kyr. E'tot is the 
total lepton energy between E'min and E'max- 





V 


P 


TP 

-•-^min 


E 

-"-^max 


'^max 


-E'tot 










(eV) 


(Hz) 


(erg) 


radio 


0.35 


1.78 




2 X 10^° 


lO^o 


3 X 10^8 


radio 


0.45 


1.78 




6 X 10^° 


10^1 


4 X lO^s 


X-ray 


0.003 


2.0 




3 X 10^^ 


2 X 10^0 


2 X 10^6 



T 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 r 




8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 

LOG^o(FREQUENCY/Hz) 



Fig. 1. — Dashed lines: The radio synchrotron and high energy (HE) IC 7-ray spectra 
of the total 2° x 3° Vela X structure corresponding to synchrotron cutoffs of 10^° Hz {thick 
dashed line) and 10^^ Hz {thin dashed line). HE IC 7-ray features from scattering on (a) 
the CMBR (b) FIR and (c) starlight. Thick solid lines: The X-ray synchrotron and 
VHE 7-ray IC spect ra including spin evolution. Tr iangles: Radio flux mea surements from 



Alvarez et all (120011 ). Circles: H.E.S.S. data from lAharonian et al. I (120061 ). Short thick 
sohd lines: INTEGRAL ("INT") data from L. Kuiper (2007, personal communication), 
i?e]9]5oSAX "MEGS" shown a s a lo wer limit for an extraction radius of 15 ar cmin and PDS 
("B") from iMangano et aD J2005h . Arrows: EGRET Vela X ULs from Ide Jager et aL 



(1l996al ). Crosses: Fermi LAT 1-year sensitivity from EGRET ULs, scaled down by x30. 



